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Abstract 

A ^-modified version of the central limit theorem due to Umarov et al. 
affirms that q-Gaussians are attractors under addition and rescaling 
of certain classes of strongly correlated random variables. The proof 
of this theorem rests on a nonlinear (/-modified Fourier transform. By 
exhibiting an invariance property we show that this Fourier transform 
does not have an inverse. As a consequence, the theorem falls short 
of achieving its stated goal. 
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1 Introduction 



Since over two decades work by Tsallis and co-workers [U EJ [3] has drawn 
attention to problems in statistical physics where Boltzmann-Gibbs statistics 
is not or not directly applicable. These problems typically involve long-range 
interactions and correlations. Tsallis initiated an approach to such systems 
that is often referred to as q-modified or nonextensive statistical mechanics. 
One of its characteristics is the occurrence of so-called g-Gaussian probability 
distributions; these are defined by 



where q is an arbitrary real number, C q is the normalization constant, and 
the domain of definition is the real a;- axis (for q > 1) or the part thereof 
where 1 + (q — l)x 2 > (for q < 1). For q — > 1 equation ( 11. ip reduces to an 
ordinary Gaussian. 

In g-modified statistical mechanics it is expected that g-Gaussians arise 
naturally when one considers sums of strongly correlated variables, at least 
for a certain class of correlations; this by analogy to ordinary Gaussians, 
that describe sums of sufficiently weakly correlated variables. Therefore, the 
observation of g-Gaussians, whether in nature or in numerical simulations, 
would lend support to the applicability of g-modified statistical mechanics. 

There have recently been attempts to provide a theoretical basis for this 
special status attributed to g-Gaussians. One attempt has consisted in nu- 
merically generating and adding up random variables with strong correla- 
tions of some well-controlled kind, and determining the distribution of the 
sum. In two examples studied by Thistleton et al. jH [5] and by Moyano et 
al. [6] the sum distribution is numerically virtually indistinguishable from 
a g-Gaussian and each was initially believed by their authors to be one. In 
subsequent analytical work Hilhorst and Schehr [7j showed, however, that 
in fact the analytic expressions of these sum distributions are unrelated to 
g-Gaussians. It was pointed out [S] , moreover, that correlations between 
variables may be engineered in such a way that the sum of N of them, scaled 
and in the limit iV — > oo, has any desired distribution. 

The question therefore is not whether a sum of strongly correlated vari- 
ables can have a g-Gaussian distribution^, but whether under the operation 
of addition and scaling of random variables such g-Gaussians appear as at- 
tractors, at least for certain classes of correlations wide enough to be of 
interest for physics. 

1 This is indeed possible; among the examples are those of references 018]. 
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A second attempt to mathematically consolidate the special role attri- 
buted to the g-Gaussian probability law is due to Umarov et al. [10J. These 
authors, in the wake of earlier conjectures [HI [6], prove a g- modified central 
limit theorem (g-CLT) in which the limit functions are q- Gaussian^. We 
recently expressed concern [5] about two points in the proof of this theorem: 

(i) First, let x\,x 2 , ■ ■ ■ ,x n , . . . denote the sequence of random variables 
whose scaled partial sums we wish to study Then in order for the g-CLT 
to apply, the first iV of these variables, for all N = 1,2,..., must be cor- 
related according to a certain condition stated explicitly in [10] and termed 
u q- independence" . This property indeed holds, according to plj] , for a conve- 
niently chosen sequence of random variables that are themselves, individually, 
g-Gaussian distributed. However, to our knowledge, no other example that 
would illustrate the theorem has been exhibited so far. It is therefore not 
sure that it is at all possible to fulfill the conditions of the theorem in any 
nontrivial way. 

(ii) Secondly, the proof of the theorem is based on the use of a (/-modifi- 
cation of the ordinary Fourier transform. This g-FT is nonlinear and, again, 
tailored to g-Gaussian distributions: when applied to a g-Gaussian, it pro- 
duces (up to scaling factors) a g'-Gaussian where g'(g) is a known function. 
It was observed in reference [Sj that this g-FT does not have an inverse. This 
observation is the subject of the present note. 

We will first show here the proof of noninvertibility hinted at in reference 
[8] and which is based on a counterexample. Whereas a single counterexample 
may be mathematically sufficient to make a point, we next ask whether one 
can perhaps obtain an invertible transformation by restricting the domain of 
action of the g-FT to a suitable subspace of probability distributions. The 
answer turns out to be negative. By explicit construction we exhibit families 
of functions all having the same g-Fourier transform; and we show that the 
g-Gaussians themselves are part of such families. 

Sections [2] and [3] deal with the definition and the non-invertibility of the 
g- Fourier transform, respectively. In section H] we present a brief discussion 
and conclusion. 

2 A (/-modified Fourier transform 

Let f(x) > be an integrable function on the real axis. Umarov et al. [TO] 
define its nonlinear q- Fourier transform / g (£) as 

= f°°dx ^ -3- , (2.1) 

J-oo [l _ ( q - tyifrfi-^x)]*- 1 

2 Extensions of the theorem to the multivariate case [T2] and to q-modified a-stable 
Levy distributions [T3] have appeared since. 
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where £ is real and from now on q > 1; furthermore here and henceforth 



./ 



x 



[f{x)] q 1 . For g — t- 1 + expression ( 12. ip reduces to the standard 



Fourier transform. If f(x) is a probability distribution, as will be the case, 
the normalization condition 



dx f(x) 



1 



(2.2) 



has to be imposed. 

The g-FT (12. ip has since its introduction been discussed on various occa- 
sions by the original authors and others [T3J [T5J [TBI [T7J [TS] . In what follows 
we will show that it is not invertible, even when restricted to the space of 
probability distributions. 



3 Noninvertibility of the q-FT 

3.1 An example 

We begin by examining a particular case ([8], footnote [30]). 



Example 1. Let us consider 
/(*) = 



xJ 




x E (a, b), 
otherwise, 



(3.1) 



where < a < b and A > 0. Due to the normalization (12. 2ft the constant A 
can be expressed in terms of the interval end points a and b, which yields 



A 



q-2 



(6 9-1 



a 9-1 J 



-(9-1) 



(3.2) 



Substituting ( 13. ip in (12. ip and doing the integral using the normalization of 
/ leads to 

= -t: • (3-3) 

[l-(g-l)i£A] 

This is the same result for any interval (a, b) that satisfies ( 13. 2 p with a fixed 
A. Hence equations ( I3.ip - (l3.2p define a one-parameter family of normalized 
functions f(x) all having the same q- Fourier transform ( 13. 3p . This coun- 
terexample shows that the q-FT is not invertible on the space of probability 
distributions. 



4 



3.2 An invar iance property 

For physicists it is important to know that the above example is not an 
isolated case that might somehow be eliminated by suitably restricting the 
space of functions. We will therefore show how other examples may be con- 
structed, many of them corresponding to physically reasonable probability 
distributions. The construction rests on the fact that for a large class of 
functions f(x) the expression 

X( x ) = xf q -\x) (3.4) 

is not invertible to a single- valued function x(X). Let us rewrite the q-FT of 
equation (12. ip as 

UO = dA F(A) -3- (3.5) 

J-oo [1 — — 1) i^A] <J - 1 

in which 

/oo 
dx 5(xf q -\x) - X) f(x) 
-Or 

= £ 

i 

where the sum runs through the set {xi} of the solutions of xf q ~ 1 {x) = X. 
Obviously, if two distinct functions fi(x) and f2(x), when substituted in 
(13.61) . lead to the same F(X), then by ( 13. 5ft they will have the same q-FT, 
This invariance leads to the nonuniqueness of the inverse q-FT. 

3.3 A class of symmetric functions 

An investigation of this invariance in its most general form would probably 
begin by classifying the f(x) according to the number of terms that they 
engender in the sum on i in equation (I3.6p . It is not needed for our purpose 
to embark on so broad an enterprise. We will study here a more limited 
but important class of functions f(x) which contains, in particular, the q- 
Gaussians. None of the restrictions on f(x) to be adopted below is essential; 
each can be overcome by a little more work. 

Let us consider for convenience a symmetric function, f(x) = f(—x), 
that is finite in the origin, /(0) < oo. Its values for x > (for x < 0) then 
determine F(X) for A > (for A < 0) and we have F(X) = F(—X). We 
may therefore limit our analysis to x, X > 0. Let us furthermore restrict our 
attention to those f(x) for which X(x) = xf q ~ l (x) is monotonic on (0, oo) 



(3.6) 
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except for passing through a single maximum at x = x m ; and let A m = X(x m ). 
We note that the monotonicity condition is not a very severe one and is 
satisfied, in particular, by the g-Gaussians (11. ip . 

Let now f(x) be in the class delimited above. Full generality not being 
our purpose, let us suppose that 1 < q < 2. In that case the integrability of 
f(x) implies that we must have 

A(0) = A(oo) = 0. (3.7) 

Let us denote by x_(A) and x + (X) the inverses of X(x) on the intervals [0, x m ] 
and [x m , oo), respectively. The pair x±(X) is an alternative representation of 
f(x). Using that dx + /d\ < < dx^/dX and that f(x) = (X/x)~ we obtain 
from (I3.6p for F(X) the expression 

which we cast in the final form 

F(X) = fEf [xr 1 (A) - xf 2 (A)] (3.9) 

where < A < A m . Starting from ( 12. 2p and replacing x by A as the variable of 
integration, we find along similar lines that the normalization of a symmetric 
function f(x) can be expressed as 

/■Am 

2 / dAF(A) = 1. (3.10) 
Jo 

Equation (13. 9p shows that the invariance can now be expressed as follows. 
The function -F(A) does not change if in (13. 9 p we substitute x±(X) H- x±(A) 
with 

xr^(A) = x!^(A) + G(A) + GL, 

^(A) = x^(X) + G(X) + G + , (3.11) 

where G + and G_ are constants and G(X) is an arbitrary function. Equation 
(13.51) shows that this substitution does not change the g-Fourier transform 
f q (£,) and ( 13.101) shows that it does not change the normalization of f(x). The 
function G(X) should satisfy certain rather mild constraints coming from 
the fact that we want the pair x±(X) to be again the representation of a 
probability distribution f(x). This distribution f(x), if it exists, can therefore 
be constructed in the following two steps: 

1. Invert x±(X) to a function X(x). 
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2. Solve f(x) from xf q l (x) = X(x), which gives 



\(x) 



x 



/(*)= ( — P- (3.12) 



The result is a function f(x) different from f(x) but which has the same 
q-FT. 

3.4 Three more examples 

We will exhibit below three examples that result from an implementation of 
the procedure of section 13.31 We have not sought to exploit the full freedom 
offered by the occurrence of the arbitrary function G(X) in (13. lip , but re- 
placed it with a single parameter. 

Example 2. Let us take for f(x) the g-Gaussian G q (x) of equation (11. ip . 
For this function direct calculation gives 

C*- 1 ± \C 2 q [q - l) - A(q - 1)X 2 V 2 
x±(A) = -2 1 ' _ ^ } —*-, 0<A<A m , (3.13) 



with 



\ m = U q -\)-\ci l . (3.14) 



2 VI ^1 

Expression (13.131) has the properties 

{oo _i 

^ and x±(A m ) = x m = (q- 1) 2. (3.15) 

Let us choose G(X) = A where A > is a parameter, and G + = G- = in 
(13TTT]) . Then 

x ± (A) = (xr(A)+i)^. (3.16) 

In expression (13. 161) the x±(X) are given by (13. 131) and therefore the x±(A) 
can be inverted to a function A(x). The result is 



A ( £ ) = Tf^' V = (?-l) a ^-^J ■ (3-17) 

We will now indicate the A dependence explicitly and write Ja{x) for the 
function j{x) represented by (13. 1 7D . Using (13. 121) we find from (13.171) 



Ux) < X{xhq - 1 



X 



Xi 



~ [ 1 + (q - 1) (x^ - A)"- 2 ]"- 1 



q-2 - q-2 

valid in the domain A < xi- 1 < oo, whereas /^(x) = for < xi- 1 < A 
and jfU(— x) = /a(^)- We have not indicated the g-dependence of this family 
explicitly. The A = member of family (I3.18P is the original g-Gaussian 
G q (x). By construction all Ja{x) have the same q-FT, independently of A. 
Hence we have constructed a one-parameter family containing a g-Gaussian 
and for which the g-FT has no inverse. 

Example 3. We consider the special case q — § of equation (I3.18p . The 
reason is that, whereas different statements concerning the g-FT occurring in 
the literature may have different domains of validity on the q axis, virtually 
all of them apply for 1 < q < 2. Hence a value in the middle of this interval 
is among the most relevant ones. For q = | equation ( 13. 181) simplifies to 

Cz (xr 1 - A)- 2 
f A (x) = 2 . (3.i9) 

X 2 [l + i(x-!-A)- 2 ] 2 

We will return to this example in our conclusion. 

Example 4. We consider the special case q = 2 of equation f 1 3 . 1 8 1) . The 
limit q — > 2 is singular. To take this limit we replace A by a parameter a 
defined as 

A = log a. (3.20) 

q - 1 

The conditions A > and 1 < q < 2 that we imposed above now require 
that a > 1. Substituting f)3.20p for A in ( I3.18P and taking the limit q — > 2~ 
we get, with an obvious change of notation, 

/«(*) = T^1T> ( 3 -21) 
1 + a z x z 

which is equal to aG^ax). Hence we have found here that the 2-FT of 
aG2(ax) is independent of a, for a > 1. This independence is in fact valid 
for all a > and is of course easily demonstrated by direct calculation. 



4 Conclusion 

The authors of reference [10] and later work on the g-FT have certainly been 
aware of the necessity for that transformation to have an inverse. Reference 
[T4] deals exclusively with this issue. However, the authors limit their inves- 
tigation to the question of whether a q- Gaussian obtained under the q-FT 
has a unique preimage in the subspace of q-Gaussians, that is, to asking if 
g'(g) has an inverse. The answer to the question thus narrowed down is affir- 
mative but has little bearing on the problem of the invertibility of the g-FT 
on a full space of functions. 
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Our example 3 appears^ in the overview by Tsallis [18] , who pays ample 
attention there to the problem of invertibility. Tsallis proposes to select a 
specific value of the parameter A, hence a specific member of the family 
/a(^), by prescribing the second moment of that function. Our view, quite 
apart from other questions that such a procedure may raise, is that this is an 
ad hoc fix covering the one-parameter case. It does not lift the degeneracy 
in cases - whose existence we have made plausible - where a single q- Fourier 
transform is associated with a many-parameter family or with a family that 
depends on an arbitrary function, that is, on a continuum of parameters. 
The presentation of our derivation in section [3] had, precisely, the purpose of 
showing the extent of the invertibility problem. 

In summary, we have shown that, when considered on a reasonably large 
space of functions, the g-modified Fourier transformation employed in the 
work by Umarov et al. JTU] does not possess an inverse. As a consequence, 
that work remains unsuccessful in its attempt of showing that g-Gaussians 
are attractors under addition and rescaling. 

The numerical search for g-Gaussians has dealt with statistical models 
that are less amenable to analytic treatment than those of references [U El 
EJ [7] mentioned above. This search has in particular focused on the logistic 
map and the Hamiltonian Mean Field Model (see [8] for references), where 
various types of scaled partial sums have been proposed as candidates for 
g-Gaussian distributed variables. In the future still other models will no 
doubt be examined with the same purpose. We leave numerical work of that 
nature out of the present discussion because, first, strictly mathematically 
speaking we have nothing to say about it; and, secondly, the interpretation of 
the numerical results has in each case been controversial. The only relevant 
corollary of the present note is that fitting numerical data with g-Gaussians 
cannot be justified on the basis of a g-central limit theorem. 
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